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Abstract 

The properties of hard photon radiation in W pair production at LEP 2 are stud- 
ied, with emphasis on the energy loss relevant to the W mass measurement. We 
use a combination of the exact 1-photon matrix element and leading logarithmic 
structure functions. Defining unobservable, observable and initial state photons in 
the phase space, it is shown that neither the one-photon matrix element nor the 
leading logarithmic structure functions alone give an adequate description of the 
energy loss due to observable or initial state photons. An event generator based on 
these calculations is available. 
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1 Introduction 



One of the major goals of the LEP 2 experiments is an accurate determina- 
tion of the mass of the W boson. A comparison of the measured value with 
the one predicted from other precision electro-weak measurements will serve 
as a sensitive check of the Standard Model. The most precise measurement 
method seems to be a constrained fit of the final state in the semi-leptonic 
(lepton+jets) and hadronic (4 jets) channels. These fits assume that there is 
no hard radiation or that the photon spectrum is known. Final state radia- 
tion can easily be incorporated by adding the photon to the nearest outgoing 
charged particle, but the shift in the fitted W mass caused by initial state 
radiation will have to be determined using theoretical input. In particular, in 
order to get a precision in M w better than 10 MeV, the average energy loss 
has to be known to 20-25 MeV. 

In W pair production the distinction between initial and final state radiation 
is not unique (unlike in the case of neutral particles such as the Z boson). 
In the matrix element, the universal leading logarithmic parts (proportional 
to \og(s/m 2 e )) are easily separable, but the non-universal finite terms do not 
split naturally (although the current splitting technique [1] seems to achieve 
a 'good' separation). We therefore define initial state radiation in the phase 
space rather than the matrix element; this is also more in line with the ex- 
perimental possibilities. We use two definitions. The first is the criterion that 
the outgoing photon is closer to the beam than to any charged final state par- 
ticle. The second assumes that events with observable photons are analysed 
separately, and defines initial state radiation as unobservable photons closer 
to the beam than to other charges particles; this is almost equivalent to a cone 
around the beam. 

The first computation of hard radiation to off-shell W pair production was 
performed by Aeppli & Wyler [2]; this calculation did not include mass effects 
so it could not cover the collinear regions. In Ref. [3] we gave the extension 
to the full phase space and an event generator. Other event generators con- 
tain a leading log description of initial state radiation, like [4-6], where the 
last two use a shower algorithm to generate px and multiple photons. The 
semi-analytical calculation of Ref. [1] also gives some of the finite terms in 
the current splitting scheme, which allows for a gauge-invariant separation 
of initial and final state separation. They include the universal non-resonant 
graphs. Finally, a first exploration of the hard radiation off the single W (t- 
channel) graphs of e + e~ — > e~u e ud^ was presented in Ref. [7]. The goal of 
this article is to give a very good description of hard radiation, including the 
exact I7 matrix element at large angles, and the leading 0(a 2 ) contributions 
at small angles to the beam. The various definitions of initial state radiation 
are compared, and the validity of the approximations made by other compu- 
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tations investigated. The influence of the universal non-resonant diagrams is 
confirmed to be small; other diagrams leading to some of the 4-fermion + 7 
final state, such as ZZ graphs, single W production and QCD graphs, have 
not yet been included. 

We start with a technical description of the matrix elements and phase space 
used. Next we give results for the energy and pt spectrum of observable, ini- 
tial state and unobservable radiation, and the average energy loss. This is 
compared with the same quantities in the one-photon and leading log approx- 
imation, and the influence of the universal non-resonant graphs is computed. 
Finally, we mention where the event generator which incorporates these cal- 
culations can be obtained. 



2 Method 



2.1 Hard Bremsstrahlung 



The method used in this calculation is an exact matrix element for the reac- 
tion e + e~ — > 4f + 7 convoluted with structure functions to describe multiple 
collinear photons. The phase space is thus the one-photon phase space (the 
multi-channel MC described in Ref. [3]) plus collinear photons, the sum of 
which is represented by one particle in each direction. No attempt has been 
made to use a shower algorithm to resolve these into individual photons. 

We use two matrix elements. The first one is the original massless one of 
Ref. [2] (explicit expressions for this matrix element are also given in Ref. 
[8]). Leading mass effects are added as described in Ref. [3]: When the angle 
between the photon and a charged particle % is so small that the mass effects 
become significant we use 

\M I 2 — \M 12 (QiQ-r) 2 m i 

1 1 vl massive | | 1 vl massless \ , \ ( \2 \ ' 



where the qj denote the massless momenta, and the pj the massive ones. This 
way the logarithmically mass singular and finite (mf/mf) terms are included 
correctly. 

The second matrix element, used mainly for cross-checks, has been generated 
automatically using MadGraph [9]. It is an order of magnitude slower but 
contains all mass effects without approximations. Using this matrix element 
one also has the possibility to study the effect of universal non-resonant graphs. 
The ZZ and single W graphs are also included, but as long as the appropriate 
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channels have not been added to the phase space mappings these can not be 
used.Q 

Given the number of expected events these I7 matrix elements are accurate 
enough for a good description of large-angle hard photon emission. However, 
they fall short in describing the collinear initial state bremsstrahlung, which 
is dominated by large logarithms log(s/m.g). These can be resummed (and 
higher order contributions can be included) in structure functions (transition 
functions) [10-14], but only after integrating out the transverse momentum 
Pt- As in Ref. [15], we combine these two approaches by only resumming the 
radiation inside a cone 9 < 6 C ; subtracting the leading log 0(a) part from 
the explicit matrix element to avoid double counting. This gives multi-photon 
emission, but only on the side where the explicit photon happened to be. 

Using structure functions one would like to write the matrix element as 



with fc ± the fraction of the incoming momentum taken away by collinear initial 
state radiation off the e ± and M(s) some hard scattering amplitude. We use 
the structure functions given in Refs [1O,16]0 

p(k) = (3k^(l + + 5^ +s )+5f(k) + 8*(k) , (3) 



with (3 = —{L— 1) and L = log(/i 2 /m 2 e ) the large collinear logarithm. The terms 
S denote the infrared finite parts of the leading one- and two-loop corrections: 



a ,3 



7T 



6^ = —(~L+ — -2) 
1 2tt V 2 3 ; 

5?(k) = ^-(l-L)(2-k) 



xV+S 
°2 



*?(*) = (f 



a a/9 ffV 2 

1 + (1-A;) 2 



log(l - k) 



+ (2-fc)(ilog(l-fc)-21og(*)-|; 



k 



(4) 
(5) 

(6) 
(7) 



2 A possible exception is the study of large-angle electrons in the case of single W 
production, as in Ref. [7]. 

3 The difference with the 3-loop YFS structure function [17] is 0(1O -3 ) in the 
variables studied. 
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For the scale we use p? = s(l — cos# c )/2, which reproduces the 0(a) results 
and agrees with the scale usually taken when no cone is defined. 

We have to define the hard scattering amplitude in such a way that collinear 
photon radiation is not included twice, preferably keeping it positive definite. 
The procedure followed resembles the one described in Ref. [15]. For large 
angles (larger than 6 C ) we take for this amplitude the I7 matrix element. 
However, inside the cone we must first subtract from this the 0(a) leading 
log part which has already been included in the structure functions: 



Here k = E 1 /E e and |Mo| 2 denotes the matrix element without photon eval- 
uated at a lower energy. This expression only holds strictly in the collinear 
limit, and breaks down when the cone is chosen too large (typically around 
90°). 

Next we add again the exponentiated leading log result: the lowest order ma- 
trix element times structure function, but now at k' = k ± + x, with x = E^/E e 
the fraction of the beam energy taken by the explicit photon and the direction 
(±) determined by the direction of the explicit photon: 



|M add ((l-A; + )(l-Ar) S )r^ c 

I |M(°>((l-ArF)(l-fc'))| 2 p(k') g(x,k>) 
J 1-k' p(k±) Jl ;JE Jx'g(x',k')- 



^min/ 



g(x, k') is an arbitrary function and S a symmetry factor defined later. In 
order to get roughly the same behaviour as the original matrix element we 
choose 

g(x, k') » Y^T' p{k ' ~ x) P2(;r/(1 ~ k ' + x)) ' (10) 



with p2{x) the function p(x, a — 2a) which approximates the matrix element 
squared in one direction. We use an approximation as the function has to be 
analytically integrable. 

However, this is only half the story, as the lower bound on the photon energy 
is now always imposed on the side on which the explicit photon happened to 
be. The region in phase space where the collinear photon on the other side has 
satisfied this bound, k T E e > E min , is included by defining S = 2 — 0(E min — 
k T E e ). The step function avoids double counting when both collinear photons 
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have large energies. The total matrix element is then given by 

|M full | 2 = |M massive | 2 - |M sub | 2 + |M add | 2 (11) 

The collinear energy lost in the forward and backward directions is the sum of 
many photons. We can approximate the px of this bunch of collinear photons 
by the pr spectrum of a single photon of the same energy to obtain a smooth 
transition from the I7 matrix element to the exponentiated one. The resulting 
violations of energy-momentum conservation are of order E~, sin 2 9, and hence 
negligible for small cone sizes. The implementation of a true parton shower 
algorithm [5] is foreseen for the near future. 

The same subtraction of the leading logarithmic contribution can be performed 
for final state photons. As the exact difference between a single fermion and a 
fermion accompanied by collinear photons is not important for most measure- 
ments, we have not performed the exponentiation as we did for the initial state 
photons. Also, radiation off a hadronic final state is badly described by the 
radiation of an on-shell quark: the effective mass of the radiating particle is the 
mass of the whole jet, not the on-shell mass of a naked quark. This reduces the 
bremsstrahlung considerably. In the time domain this means that the quark 
will radiate off gluons before it can radiate off a photon; the amount of hard 
radiation near a jet is thus overestimated. The possibility to turn off the final 
state radiation means that one can leave the generation of these photons to 
more heuristic algorithms (see, e.g., [6]) tuned to the observed spectrum if this 
spectrum is important. 

An unsolved problem is the matter of gauge invariance of bremsstrahlung off 
unstable particles. As we use a constant width in our calculations the current 
is conserved (the substitution p^ — > gives zero), but this is not sufficient to 
guarantee gauge invariance. However, there are no cancellations which could 
amplify the gauge violating terms, so we expect the effects to be small. A test 
using the inclusion of triangle graphs is being performed. 

2.2 Soft Bremsstrahlung 

In the following we will consider mainly the average energy lost by photons, 

{ B,) = £fdE,E,£. (12) 

The normalisation o"q r is computed from the tree level matrix element for the 
reaction e + e~ — > 4 fermions, convoluted with structure functions. This has an 
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accuracy of a few percent. However, as an overall factor it will not influence 
the conclusions very much. 

In the case of a strict one-photon calculation the integral in Eq. 12 is not sen- 
sitive to soft photons, and we can make accurate predictions without knowing 
the virtual and soft corrections. However, when we add the possibility of mul- 
tiple photon emission through the procedure given above, and define E 1 as 
the sum of the photon energies, these do play a role. Note that we consider 
the regime E min <C T w . In this case the soft corrections are easily found to be 

|M soft | 2 = e 2 ^2Q i Q j B{ Pi ,p j ;E Ilda ,X)\M \ 2 , (13) 

i,3 



with Qi the charges of the particles in terms of the unit charge e. The function 

-C/7<-C/ m in 



is given in Ref. [18] (the sign depends on whether the particles are incoming 
or outgoing). It depends logarithmically on the infra-red regulator, for which 
we use a small photon mass A. We can include the corresponding infra-red 
logarithms of the virtual corrections by the substitution A — > /iiR. We choose 
for this scale /iiR = Mw/n, with n the number of charged final state particles. 
With this heuristic, the difference in the total cross section between the full 
result and the leading logarithmic approximation is less than 1.5% in the 
region 175 < y/s < 205. 

Both the soft and virtual corrections will contain logarithms in the virtuality of 
the W bosons, \og(p 2 — M w +iM w T w ). However, these cancel in the sum. This 
can be seen in the hard photon spectrum (Figs 1-3) as an absence of a kink at 
E 1 m Tw The cut-off on the initial-final and final-final state interference as 
the energy increase is thus exactly compensated by the emergence of radiation 
off the W bosons. 

We also extract the Sudakov double logs 

|M sud | 2 = -^E^log 2 — |Mo| 2 . (15) 

When we introduce the structure functions we have to subtract again the part 
which would be double-counted: 

|M sub , soft | 2 = ^[2(L - 1) log + \h + £ - 2] (16) 
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The tree term has already been added to the hard radiation for E 1 > E m \ n , 
we need to add it only to the soft terms as 

|M trcc | 2 = (1 - 6{{k mhi - k + ){k mhi - k-))\M \ 2 (17) 

Adding all up the estimate for the non-radiative part is 

I M full)Soft | 2 = | M soft | 2 - | M sud | 2 - | M sub , soft | 2 + | M tree | 2 (18) 

The other terms in the virtual corrections do not contain large logarithms, and 
are thus expected to be of order a/n < 1%. As they only enter as corrections 
to the main effects of hard radiation in the quantities studied in this paper this 
is negligible. For the factorizable graphs, for which the off-shell result is under 
study, this reproduces the full result to a few percent in the cross section. 
A change in the scale /ii R of a factor 2 gives a difference of 3% in the cross 
section, but only 1% in the observable energy (where the effect of the virtual 
and soft terms is much smaller). We verified that the numbers obtained with 
this soft matrix eleemnt did not any more depend on the cut-off E min . 



3 Results 

We now turn to the predictions we can make for the hard photon spectrum in 
W pair production. We first define the 'LEP 2 canonical cuts' that are used to 
compute observables that can easily be compared to other calculations. These 
cuts also give an impression what an actual experiment could observe. Next 
we show the energy and pt spectrum for observable and initial state photons. 
This can be compared with the one-photon approximation on the one hand, 
and the leading log approximation on the other. Finally we discuss the effect 
of the non-resonant diagrams. 

The canonical cuts of the ADLO/TH detector [19] are given in Table 1. r's 
are not considered. The cuts on the photon separations should be interpreted 
as follows: first each photon is combined with the nearest charged particle if 
the angle to it is smaller than 5° (outgoing) or 1° (beam); combining with 
the beam means dropping it. Next all the other cuts are applied. This way 
spurious collinear divergences near cuts are avoided. 

Observable photons are defined as those which pass these cuts. Unobservable 
initial state photons are the ones cut out by the energy cut or the minimum 
angle to the beam; the sum of their energies and px can be deduced from 
momentum conservation. Finally, we define initial state photons to be the 
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5° 


5° 


5 GeV 


jet 


3 GeV 


0° 



Table 1 

Canonical cuts 

unobservable ones plus observable photons which are closer to the beam than 
to any other charged particles. 

In the subsequent part of this work we use the following parameters. We work 
in the G M scheme, but with the coupling of the extra photon determined by 
a(n 2 = 0). All masses are taken from the Particle Data Group 1994 edition 
[20]. The quark masses were taken m u = 63 MeV, = 83 MeV, m s = 215 
MeV and m c = 1.5 GeV for compatibility with the virtual off-shell calculation. 
The W mass is 80.23 GeV. The W width is computed in leading order, with 
the strong coupling constant a s = 0.117; this gives Y w = 2.081 GeV. To 
conserve the electromagnetic current we use a fixed W width; the pole position 
is computed as fj, w = — — 5/27rTy V /M w — iM w T w — iTfy/pi + i(l — 
1 / pi^Tyy I Mw ■ The Z width is neglected, and the energy of the outgoing 
photon restricted to less than E — (Mz + 5Tz) 2 /4:E to avoid the Z peak.0 
Unless otherwise noted we assumed a beam energy of 87.5 GeV. We used a 
cone angle of 10°. Changing this to 5° shifts the observables by an amount 
comparable to the integration accuracy, 0(1%). 

The results are presented in Figs 1-3 and Tables 2 and 3. First we give the 
observable tree level non-radiative cross section (cr ), the same convoluted with 
leading log structure functions (<Jq t ) and with non- resonant graphs (o"o+ nr ). 
For the next entries we consider the full calculation described in section 2, 
the exact one-photon matrix element (I7) and the leading log result (LL). For 
technical reasons the leading logarithmic approximation can only be computed 
for photons emitted in the forward direction (6 C = 90°), so the scale is taken 
to be fi 2 = s/2. This will underestimate the leading log result. We also give 
the full result including the universal non-resonant diagrams (full+nonres). 
The statistical errors on the cross sections are O(0.1%), but the differences 
(LL — full) and ((res + nonres) — res) were computed directly in this form and 
have relative errors of a few percent and a few tens of percents respectively 
on these differences. The statistical errors on the the average photon energies 
are slightly larger, 0.3-0.5%. We introduced an upper cut-off on the photon 
energy to avoid the Z peak, this influences only the observable average energy. 



The cross section above this cut is entirely negligible [2]. 
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For the normalisation we need the total cross section for W pair production. 
As long as the full off-shell one-loop result is not yet available we use the 
lowest order cross section with leading logarithmic structure functions. As the 
final state corrections cancel against the corrections to the width, and the 
initial-final interference is expected to be of order aTw this is a reasonable 
approximation. This uncertainty only affects the overall normalisation of the 
results, not the differences between the different results. The main effect of the 
initial state radiation is of course to lower the result, as we are just above the 
threshold for W pair production and phase space is limited. The contribution 
from the non-resonant graphs is negligible in this energy range (as had already 
been observed in Ref. [2]). With the scale choice used, the full result deviates 
from this estimate by less than 1.5%. 

One can see from the Tables that an appreciable fraction (around one quarter) 
of the events will be accompanied by photons observable in the ADLO/TH set 
of cuts. This is due to the excellent forward coverage (0 7 > 1°) and electromag- 
netic calorimeter (E^ > 0.1 GeV) assumed in the canonical cuts. Reducing the 
angle to 10° this fraction still is around 20%, of which half also has an E 1 > 1 
GeV. For the W mass measurement it would be advantageous to make use 
of the extra information present and analyse these events separately. Neither 
the I7 matrix element nor the leading log approximation give a satisfactory 
description of the observable photons. By radiating only one photon one ef- 
fectively normalises to the non-radiative lowest order, which is too large. On 
the other hand, the leading log approximation misses the negative initial-final 
state interference terms {E 1 < rV) and the radiation off the (off-shell) W 
bosons (mainly when E 7 > rv).Q Finally, given that the un-exponentiated 
large-angle contribution of the cross section still is 20% of the total cross sec- 
tion the region E 1 < 1 GeV should not be trusted to 1% even in the full 
calculation. Most of the cross section here is, however, associated with the 
final state and hence does not influence the W mass measurement. 

The observable photon energy is dominated by final state radiation and hence 
not very interesting. The unobservable energy spectrum is much more indepen- 
dent of the final state. It is not completely independent due to the possibility 
of observable jets in the beam pipe (8j = 0°). The initial state radiation as- 
sociated with these events is sometimes associated with the final state by the 
canonical cuts, thus lowering the average energy. As the radiation off jets is 
not modelled correctly anyway, a cut will have to be imposed to exclude this 
contamination. One sees that analysing observable photons separately reduces 
the average energy loss, and hence the size of the theoretical corrections to be 
applied to the fitted W mass. The difference in (E™ 1 ) is due to the amount of 



5 The fact that there is no discernible structure at £L ~ Tw in Fig. 1 suggests 
that these two contributions are intimately connected; indeed, one cannot define 
the 'radiation off the W in a gauge invariant way. 
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phase space available: 2, 3 or 4 charged particles. For the unobservable radia- 
tion the leading log approximation is, as expected, quite goodQ. As in the total 
cross section, the contribution of the non-resonant graphs is not visible due 
to the limited integration accuracy. The contributions from the universal non- 
resonant graphs cannot be distinguished from the statistical fluctuations,^ 
but may just be relevant. At 205 GeV the integrated radiative corrections are 
smaller, but as there is much more phase space for the photons, more can be 
observed. The same remarks on the accuracy of the various approximations 
hold here, except that the deviations tend to be larger. 



4 Conclusions 



We have studied the properties of hard radiation in W pair production at LEP 
with a combination of the exact matrix element for one-photon emission and 
resummed structure functions for the initial state collinear large logarithms. 
We find that a large fraction of the W events will be accompanied by ob- 
servable hard photons; a separate analysis of these events should decrease the 
systematic errors in the W mass measurement. These events are not well de- 
scribed either by the I7 matrix element or a pure leading log approximation. 
The treatment of photons near jets is an unsolved problem. 

The average energy lost by the remaining, unobservable photons is smaller. 
The canonical LEP 2 cuts leave a contamination of final state radiation due to 
the inclusion of jets down to the beam pipe, which will have to be removed. In 
this region of phase space, a leading log treatment is adequate. We have not 
yet addressed the issues of radiation off non-resonant graphs other than the 
(probably negligible) universal ones; in particular the QCD, ZZ and single W 
graphs have not yet been included. 

An event generator based on the calculations presented in this article can be 
obtained from [ftp : // rulgm4 . LeidenUniv . nl/pub/gj| or via the World Wide 
Web. 



The systematic shift for leptonic, semi-leptonic and hadronic channels is probably 
due to the choice of the scale //ir, which should be slightly different for the different 
final states considered. 

The numbers presented here are based on 2.5 10 5 weighted events per channel 
and took about 100 hours in all to compute on a fast workstation. 
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E heam = 102.5 GeV 





leptonic semi-leptonic hadronic 


cr [pb] 
<[pb] 

4+nr[pb] 


.831 (7%) 5.280 (8%) 8.316 (8%) 

.774 4.908 7.724 

.774 (.07%) 4.912 (.08%) 7.730 (.08%) 


O-obs^isr 

full 

1 7 

LL 

with nonres 


.295 .289 .267 
.354 (20%) .326 (13%) .302 (13%) 
.317 (7.2%) .303 (4.9%) .276 (3.3%) 
.296 (.08%) .289 (.08%) .268 (.02%) 


(E° hs ) [GeV] 
full 

17 
LL 

with nonres 


2.010 1.802 1.555 
2.106 (.096) 1.830 (.028) 1.593 (.038) 
2.084 (.074) 1.853 (.051) 1.599 (.044) 
2.008 1.801 1.554 


(£Junobe) [ Ge y] 

full 

17 

LL 

with nonres 


1.902 1.923 1.884 
2.266 (.364) 2.240 (.316) 2.241 (.357) 
1.866 (-.035) 1.880 (-.044) 1.875 (-.009) 
1.900 1.922 1.882 


(E™) [GeV] 
full 

17 
LL 


3.168 3.147 2.976 
3.549 (.380) 3.447 (.300) 3.344 (.368) 
3.156 (-.013) 3.106 (-.041) 2.992 (.017) 



Table 3 

Cross section for observable photons and energy lost to observable, unobservable 
and initial state photons at y/s = 205 GeV. 



14 



Figures 



Inclusion of the universal non-resonant graphs gives curves which are indis- 
tinguishable from the 'full' calculations. The wiggles are due to the limited 
statistical accuracy (10 7 weighted events per curve). 
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Fig. 1. Observable photon energy in the semi-leptonic channel. 
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Fig. 2. Unobservable photon energy in the semi-leptonic channel. 
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Fig. 3. Initial state photon energy in the semi-leptonic channel. 
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